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Abstract 

In this paper we present a new closed-form solution to a chaotic difference equation, yn+i = a2y^ + 
aiUn + ao with coefficient ao = (ai — 4)(ai + 2)/(4a2), and using this solution, show how corresponding 
exact roots to a special set of related polynomials of order 2^' ,p G N with two independent parameters 
can be generated, for any p. 
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1 Introduction 

The difference equation yn+i = fiUn), n = 0, 1, 2 . . . where 

i-/ N 2 (ai - 4)(ai +2) 

4a2 

and ai, a2 G C, 02 7^ has the exact, general solution 

y„H = -(2cos(c.2")-^). (2) 
02 V 2 / 

We discovered this solution from the solution presented in [1] by reparametrising with two new variables as, 
for example, y„ — ccos(w2") -I- d. We then inserted this into the difference equation and equated powers of 
cos(ci;2"), obtaining that c = 2/a2 and d = — ai/2a2. We were also interested to notice that the fixed-point 
problem for (1) is also a root- finding problem for the polynomial: 

Fiy) - y - (3) 

There is a result sometimes known as Abel's impossibility theorem which states that there are no exact 
expressions for finding the roots of general polynomials of order greater than four in terms of a finite number 
of elementary operations. Therefore in general we would be forced to use an iterative root-finding method for 
this polynomial (see [2] and [3] for example). These methods are typically difference equations, for example 
the method of Newton-Raphson iteration for solving the equation F{y) — is, given a close guess yo: 



(4) 



The main requirement of this method is convergence to some fixed point. However, a problem sometimes 
arises in that the system (4) may well oscillate, or in some cases behave chaotically. Af uch progress has been 
made by numerical analysts in proving convergence and inventing better methods that are stable even for bad 
initial guess values, see for example [4]. Some beautiful methods have been devised that combine numerics 
with topology. They make it possible to find approximate solutions of a given system by a continuous 
deformation of the solutions of a related one that is exactly solvable, guaranteeing global convergence [5]. 
However, in this case, as in other special cases, we can find the roots of the polynomial (3) exactly, without 
the need for iteration. The rest of the paper shows how to do this explicitly. 



2 Periodic Orbits Are Roots of the Polynomial 

The impossibility theorem is a general one: there are rare and special cases where it does not hold, and this 
paper presents another of these special cases, using a result from chaotic dynamics. The difference system 
(1) has two fixed points. It also has a countable infinity of periodic orbits of all cycle lengths. Rearranging 
the periodic orbit equation gives equation (3) which is a polynomial equation of order 2^, the first two of 
which are: 

2 , /„ iv. , (ai-4)(ai+2) ^ 



a2y + (ai - l)y + 
and 



4a2 



iai^ - 4ai - 1 ) ?y + — ( -^ai'* - ai^ + \a\ + 2 ) = 0. 

2 / a2 ylo 2 / 

Thus, in finding the periodic points of the equation /^(t/) = y, we also find exact solutions to the polynomial 
equations (3). We shall now show the technical details involved in finding these solutions. 
Since the cosine function is bounded, so too is the solution (2). In addition, as n increases, the binary 
expansion of the expression t<;2" shifts successively leftward, as described in [1], and many other texts on 
chaotic dynamical systems (see for example [6] or [7]). A key result is that if the binary digit expansion of 
a;/27r is periodic, then so too is the behaviour of the solution. Thus, finding the iojl-K that have periodic 
binary expansions leads us to the periodic points of /^(y) = y, which in turn can be interpreted as the roots 
of the 2^" order polynomial. 

For the purposes of this paper, therefore, only values of 

(/> = a;/27r 

that have periodic binary digit expansions are relevant to us. Since irrational numbers are not periodic, 
we must choose </) rational, i.e. we want 4> = k/L with G N such that </> has a periodic binary digit 
expansion, with period p = log2(m), where m is the order of the polynomial for which we wish to find the 
roots. 



3 Constructing (j) 

Rational fractions with periodic binary digit expansions with periodicity p correspond to the following 
expressions: 

oo , 
m=l 
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As mentioned in the previous section, the periodicity of the binary expansion of <j) imphes periodicity of the 
solutions to (3). To see this, note that: 



= 2P\^ fc(2P)-'"j 

OC OO 

= (2P)fc(2P)-™ = Y A:(23')(-'"+i) 



m=l 



m=0 

oo 



m—1 



m=l 



(/)(mod 1) 



Therefore, 2P(p = 0(mod 1), which imphes, since co = 27r0, that 2^0; = a;(mod 2tt), and we reach the 
conclusion that cos (2^0;) = cos(a;), as required. 

There is, however, a minor complication to this scheme due to the symmetry of the cosine function about tt 
in the general solution (2). This symmetry implies that values of a; that lie equidistant from tt produce the 
same solution. Then 



|a; — 7r| = |(^27r — 7r| = tt 



2k 



2P - 1 



2P - 1 



|(2fc + l)-2f| 



must be unique for all choices of k. However, 



|(2n+l)-2f| = |2f-(2n+l)|, (5) 

so that k = m and fc = 2** — (m + 1) for m = 0, 1 . . . 2^ — 1 lead to symmetrically identical solutions to the 
polynomial. Therefore, only A; = 0, 1 . . . 2p~^ — 1 give the unique required solutions. 

The consequence of this is that in order to find all solutions of the 2^ order polynomial, we must seek solu- 
tions to the polynomial of order 2^p instead, since all periodic points of f^{y) = y also satisfy P'^{y) = y. 
We therefore set L — 2^^ — 1. However, the converse is not true: not all the periodic points of P'^{y) = y 
are periodic points of /^(y) — y. Therefore, in enumerating the solutions of the polynomial, (p must satisfy 
two criteria: 

(i) As mentioned above, the symmetry of the cos function implies that (j) must be less than i, which in turn 
implies that k < 2^p~^ 

(ii) (j) must either have a periodic binary digit expansion with period p, or period 2p but have uj symmetric 
about TT under the iterative shift of p digits, i.e. cf) and ^2^' lie equidistant from tt. 

Therefore, to construct and enumerate all with binary digit expansions of period p when embedded in a 
sequence of 2p, we set: 

k = m + m2P = m{2P + 1) 

for m — 0,1 . . .2P~^ — 1. We thus enumerate half of the required solutions. Secondly, from (5), to find all 
the symmetric <p with binary digit expansions of period 2p, we then choose: 

k = m2P-m = m{2P - 1) 

for TO = 1,2... 2P~^, and we have thereby enumerated the remaining solutions. 
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4 A High-Order Example 



Here we demonstrate an application of the method to the order eight polynomial P{y) — y = 0. The following 
expressions for the coefficients in ascending order of y are: 



3 1 
-1 - 16ai + lOaf - -al + —al 
2 16 

02 (-lQ + 30al-^ai + ^al 



/ 35 
all 20- 30o? + 

ai (-24oi + 7a?) 
al (-8 + 7af ) 

4aia2 



We use the LaGucrrc root-finding method [3] with a2 = — 1 + i, a\ = 2 — i, and we find that numerical 
solutions to the polynomial accurate to four decimal places are: 

y = {-0.25 -0.75i, -0.016 -0.516i, 0.1265 -0.3735Z, 0.5764 -|-0.0764i, 
0.97252 + 0.4725i, 1.25 + 0.75i, 1.651 + 1.151i, 1.6897 + 1.1897i} 

For this polynomial, p = 3 and so i = 2"^^ — 1 = 63. Next, we enumerate the solutions, firstly for 
m = 0, 1 . . . 2P-1 -1 = 3. Therefore, 



k = m(2P + l) =0,9,18,27 



or in binary notation 



9 18 271 r 1 2 3 
63' 63' 63' 63 J 1 '7' 7' 7 



= {0.000000, 0.001001, 0.010010, 0.011011} 



where for example 0.010010 indicates the infinite binary digit repetition of the sequence 010010. 

Prom these values of (f)m we then calculate the first half set of solutions, here given to eight decimal places: 

ym = — f2cos((/),„27r) - ^) 
a2 V 2 / 

= {-0.25 - 0.75i, 0.1265102 - 0.3734898i, 0.97252093+ 0.47252093i, 
1.65096887+ 1.15096887i} 

Secondly, we enumerate the symmetric set for m = 1, 2 . . . 2p~^: 

m{2P - 1) = {7, 14,21,28} 
r 7 14 21 28} r 1 2 1 41 
\63' 63' 63' 63 J \ 9' 9' 3' 9 J 
{0.000111, 0.001110, 0.010101, 0.011100} 
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from which we can calculate the second set of solutions, again to eight decimal places: 

Vm = ^ (2cos((?!)™27r) - 

= {-0.01604444 - 0.51604444i, 0.57635182 + 0.07635182i, 1.25 + 0.75i, 
1.68969262+ 1.18969262^} 

5 Conclusions 

By finding a general solution to a discrete, chaotic system and interpreting a problem of finding roots of 
a high-order polynomial as the fixed point equation for that chaotic system, we have shown how to obtain 
exact roots to the polynomial. Wc then compared this with the results of a numerical root-finding method. 
Of course we have not presented a general root-finding method, but we find it intriguing to notice that 
iterative root-finding methods are often difference equations in their own right, for which general solutions 
of particular cases may well be known. 
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